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Abstract. We present the theory of an atomic gas in an optical lattice near a 
Feshbach resonance. We derive from first principles a generalized Hubbard model, 
that incorporates all the relevant two-body physics exactly, except for the background 
atom-atom scattering. For most atoms the background interactions are negligible, but 
this is not true for ^Li, which has an exceptionally large background scattering length 
near the experimentally relevant Feshbach resonance at 834 G. Therefore, we show how 
to include background atom-atom scattering by solving the on-site two-body Feshbach 
problem exactly. We apply the obtained solution to ^Li and find that the background 
interactions indeed have a significant effect in this case. 
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1. Introduction. 

The field of ultracold atoms lias recently seen many exciting developments in which 
strong interactions play a crucial role. Examples are the observation of the quantum 
phase transition between the superfiuid and the Mott-insulator phase of a Bose gas 
in an optical lattice [T], the observation of the BEC-BCS crossover in a Fermi gas 
[21 ini m IHl El E! , and the observation of fermionization in an one-dimensional Bose gas 
[HI E] • To a large extent these developments are due to the successful implementation 
of two new experimental techniques, namely the use of optical lattices ^U] and the 
application of Feshbach resonances in the collision of two atoms fUlEI- Both techniques 
have contributed greatly to the unprecendented control over the relevant physical 
parameters of an atomic gas. 

In this article, we focus on the combination of these two important techniques and 
present the microscopic theory of an atomic gas in an optical lattice near a Feshbach 
resonance. This problem was already studied in ref. [IB], where it was argued that 
an optical lattice could be very well suited to overcome the experimental difficulties 
in studying the quantum phase transition that occurs in an atomic Bose gas near a 
Feshbach resonance ^1 . The observation of this quantum Ising transition between a 
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phase with only a molecular condensate and a phase with both an atomic and a molecular 
condensate is complicated by the fast vibrational relaxation of the Feshbach molecules 
due to collisions, which appears to prevent the creation of a molecular condensate 
In an optical lattice with low filling fractions, collisions can essentially be neglected and 
this problem is expected to be much less severe. Furthermore, in a subsequent study, 
a resonantly-interacting Bose- Fermi mixture in an optical lattice was considered, which 
led to the prediction of another XY-like quantum phase transition, associated with the 
Bose- Einstein condensation of the bosons in the mixture [T7j . 

In both studies OCZI, a generalized Hubbard model was obtained by incorporating 
the relevant two-body physics exactly. As a result, simple mean-field techniques could be 
used to accurately describe the many-body physics of the atomic gas near the Feshbach 
resonance in the optical lattice. This approach led to some discussion jl8| I19j. which 
motivated us to rederive in this article the proposed generalized Hubbard model from 
first principles by means of a field-theoretical calculation. In this manner, the validity 
of the original approach is put on an even more rigorous basis. 

Another aspect that we discuss in this article, is the effect of background atom-atom 
scattering on the Feshbach physics. This effect was not incorporated in the previously 
mentioned studies, since only systems with a small background scattering length were 
considered. However, for ^Li, which has an exceptionally large background scattering 
length near the experimentally relevant Feshbach resonance at 834 G, the neglect of the 
background interactions is no longer allowed. In this article, we therefore show how to 
include the background atom-atom scattering by exactly solving the two-body Feshbach 
problem on a single site of the optical lattice. The exact knowledge of the two-body 
physics on a single site can then be directly incorporated into a generalized Hubbard 
model that describes the many-body physics of ^Li in an optical lattice near a Feshbach 
resonance, in a similar manner as recently achieved for |20j . 

2. Effective atom-molecule theory. 

In this section, we give an ab initio field-theoretical derivation of the effective atom- 
molecule theory that describes an atomic gas in an optical lattice near a Feshbach 
resonance. The theory is formulated in the tight-binding limit and incorporates all 
the relevant two-body physics exactly. If desired, this rather technical section can be 
omitted in a first reading of the article. To facilitate this, section El has been written 
in such a way that it stands on its own. There, we show how to solve analytically 
the two-channel Feshbach problem for two atoms on a single site including background 
atom-atom scattering. We also apply the obtained theory to ^Li. 

2.1. Effective action. 

The quantity of interest in a quantum-field theory is the generating functional Z of all 
the Green's functions. This functional determines all the possible correlation functions 



Dressed molecules in an optical lattice 



3 



of the system. Specifically, let us consider the field theory for an atom-molecule gas 
that is described by the action Slipl, ip^, ip*^, ip-^^^], with ipg^ and ip^ being the atomic and 
molecular fields, respectively. The generating functional in imaginary time is defined by 

■^[■^a' "^a' "^m' '^m] 

= / dir.] rf[V'a] d[rj dbPj exp [—Sir., i^,, V'J + Sj^ (1) 

where the source currents couple to the fields according to, 

Jd^J c?x[V'*(x,r)Ja(x,r) + J*(x,r)V'a(x,r) 

r) J„,(x, r) + J^(x, r)^^(x, r)] . (2) 

By taking functional derivatives of Z with respect to the currents we can calculate all 
the correlation functions of the theory. 

Instead of working with Z we usually prefer to work with the generating functional 
W of all the connected Green's functions, which is related to Z through Z = exp{W). 
The functional derivatives of W with respect to the currents are the expectation values 
of the fields, i.e., 

SW 

= (^:(x,r))^0:(x,r) (3) 



5Ja(x, r) 
6W 

—- ^ = «Xx,r))^0;(x,r). (4) 

Similar equations hold for the expectation values of the complex conjugated fields 
'?/'a(x, r) and ■?/'m(x, r). Instead of using W, which only depends on the current sources, 
it is possible to define a functional F that depends explicitly on the fields 0a(x, r) and 
0m(x, r), and which is related to W by means of a Legendre transformation, i.e., 

0a> C 0m] = - W[J*, J^,J^, JJ 



+ 



J dr J dx [0*(x,r)Ja(x,r) + J*(x,T)0a(x,r) 
C(x, r)J,„(x, r) + j;^(x, r)(/)^(x, r)] . (5) 



The reason for defining this functional is that F is related to the exact effective action 
of the system through S'^^ = —hT. Technically, F is the generating functional of all 
one-particle irreducible vertex functions. In our case, the exact effective action for the 
atom-molecule theory can be written as 



+ ^?(CA0a + 0:0:0m) + •••], (6) 

where G*a/m is the exact propagator of the atoms/molecules, and g is the exact three- 
point vertex. The dots denote all possible other one-particle irreducible vertices, which 
turn out to be less relevant for our purposes. 

2.2. Microscopic action. 



Starting from the microscopic atom-molecule theory in an optical lattice, it is our goal 
to derive an effective quantum field theory, that contains the relevant two-body physics 
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exactly. For calculational convenience and to facilitate comparison with previous work 
[T^ . we consider the case of bosonic atoms. The generalization to fermionic atoms is 
then straightforward. 

The total microscopic action S describing resonantly-interacting atoms in an optical 
lattice can be split up in three parts, namely a purely atomic part S^, a purely molecular 
part and a coupling between atoms and molecules Sam, 

The purely atomic part Sg. is given by 

ff^p f ( h'^v'^ \ 

5'a = / rfr / dx V^*(x,r) - /i + \/o(x) Va(x,r) 

Jo J \ 2ma / 

+ / dr ^€(x, r)V^:(x, r)^a(x, r)^a(x, r) , (8) 

2 Jo J ma 

with n the chemical potential, ma the atomic mass, Obg the background scattering length 
and Vo(x) the external periodic potential due to the optical lattice. In first instance, 
we neglect the background atom-atom scattering, since we are primarily interested in 
the resonant interactions between the atoms and the molecules. For most atoms, this 
is an accurate approximation. Furthermore, the same approximation was also used in 
previous work JH] , with which we ultimately want to compare the following derivation. 
Finally, in section IHl we overcome this approximation and show how to include the effect 
of background atom-atom scattering. 

The purely molecular action Sj^ is given by 

Sra = drj C(x, r) ihdr - + SB-2fi + 2^(x) j ^,„(x, r) . 

(9) 

where Sb is the so-called bare detuning (see also sectional). The action S'am corresponding 
to the atom-molecule coupling that describes the formation of a bare molecule from two 
atoms and vice versa, is given by 

S'am = ^ dr J d:sij dx.' g{x - x') {^^((x + x')/2, r)^a(x', r) ^a(x, r) 

+V^:(x', r)^:(x, r)^„,((x + x')/2, r)} , (10) 

where g{'x) denotes the atom- molecule coupling. 

We can write the actions involving the atomic fields in a more convenient way as 

S'a + S'am = / drdr' / d:Kdx.' 

2 Jo J 

X (^:(x,r),^,(x,r))G-i(x,r;x',r') ^f^^^}^ ) , (H) 

with G the 2x2 (Nambu space) Green's function matrix, given by 

G-^ = Ga ^ - 5] = Ga-^ (1 - GaS) . (12) 



Dressed molecules in an optical lattice 



5 



Here, the atomic Green's function matrix G^^ ^ has the following form 

G'-i(x,r;x',r') 



G 



-1 







G'r^(x',r';x,r) 



(13) 



where the atomic zeroth-order Green's function Ga satisfies 



hdr 



2ma 



+ Vo{^) - G'a(x, r; x', r') = -h6{r - r')5(x - x'). (14) 



The self-energy matrix /iS is given by 





2(7(x-x')^;;((x + x')/2,r) 
X 5(r-r'). 



2(7(x-x>^((x + x')/2,r) 




(15) 



From the microscopic atom-molecule action we can calculate the grand-canonical 
partition function Z as the functional integral 



dmdmdirjdii^J exp \ -^30^1, V^,, V^n 



(16) 



Since the integral is Gaussian in the atomic fields, we can perform this integral exactly, 
giving 

Z = J d[rjd[iJm] exp Tr[ln(-G-i)] - ^sjj 

= J d[rjd[^Pj exp |-i Tr[ln(-G:i)] - i Tr[ln(l - GaS)] - ^S, 



(17) 

where the traces are taken over a 2 x 2 matrix structure, coordinate space, and imaginary 
time. By integrating out the atoms, we have obtained an effective action for the bare 
molecules that consists of three terms. The first term represents the contribution from 
the ideal atomic gas, which does not depend on the bare molecular fields. The second 
term can be brought in a more illuminating form by expanding the logarithm. Because 
of the trace, only even powers in the expansion give nonzero results, yielding 



Tr[ln(l - GaS)] 



y —(GaS)' 



n=2,4„ 



:i8) 



The first nonzero term is of second order, namely — Tr [GaSGaS] /4, and expresses 
that a molecule can break up into two atoms and recombine again. This term can be 
interpreted as corresponding to the self-energy of the bare molecules and a diagrammatic 
picture of the term is shown in figure H The higher-order terms in equation (fTH|) pertain 
to interactions between the bare molecules. 



2.3. Molecular self-energy. 

Next, we evaluate the molecular self-energy term, — Tr [GaSGaS] /4, by explicitly 
carrying out the trace. As a first step, we introduce center-of-mass coordinates 
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Figure 1. Diagrammatic representation of the term — Tr [GaSGaS] /4, which 
expresses that a bare molecule can break up into two atoms and recombine again. The 
coupling strength between the atoms and the molecule is given by g. This diagram 
can be interpreted as a self-energy of the bare molecules. 



R = (x' + x")/2, R' = (x + x'")/2 and relative coordinates r = x' - x", r' = x" 
Writing out the trace in these coordinates yields 



1 2 r^P 



/ dr dr' / dH dH' dr dr' g{r)g{r')'ip'^{Il,T) 
X V'mlR', r') Ga(R + r'/2, r; R + r/2, r')Ga(R' - r'/2, r; R - r/2, r') 
= - At dr dr' dR dRV;;(R, r)V^„,(R', r ) 

h Jo J 

X Ga(R', r; R, r')Ga(R', r; R, r'), (19) 

where in the last hne we made use of the usual pseudopotential approximation g{r) = 
g6{r). 

In order to proceed, we rewrite the atomic Green's functions by performing an 
expansion for both their spacial and temporal arguments. For the expansion of the 
space dependence, we use the complete set of atomic Bloch wave functions Xn which 
are solutions to the Schrodinger equation for a particle in a periodic potential 



+ K)(x) 



XnM = en,kX^,k(x). (20) 



2m, 

Here, the exact Bloch wave functions are given by 

X^,k(x) = -^Ee^^-'''<(^-^^)' (21) 

with tyn(x) the so-called (atomic) Wannier functions, x, the location of the optical lattice 
sites, Ng the total number of sites, k the lattice momentum, and n a set of quantum 
numbers for the various Bloch bands in the optical lattice. The Fourier expansion in 
imaginary time is performed by using the well-known Matsubara modes [21] • Combining 
the two expansions, we obtain 

Ga(x, r; x', r') = ^ E E G.{n, k; z^„)x„,k(x)x:,k(x')e-'""(^-^'^ (22) 

"■/^ n n,k 

with u!n the Matsubara frequencies and G'a(n, k; iun) given by 

-h 

G'a(n, k; iun) = — r ■ • (23) 
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Next, we also expand the space and imaginary-time arguments of the molecular fields 
?/'jn(x, r). For the temporal expansion we use again the Matsubara modes and for the 
spacial expansion we use the (molecular) Wannier functions w™(x). As a result, 

^^(x, t) = ^YY V.K)<(x - x,)e-^-"^ (24) 



with the complex conjugated version for t/'^. The operators h^^ annihilate a molecule 
with a set of quantum numbers n on site i of the lattice. 
Substituting equations and into (fTI?|l . yields 

- -Tr fG.SG.Sl ^ V- V- V- V- V- V- / 



Y.Y.Y.Y.Y.Y.j drdr'JdRdR' 

n 1-)' rv) (VT-i ' t 1 rtn m' n n' 1^' 



4 ^ 

r""^ n,n' m,m' i,j m,m' n,n' k,k' 

x&;.,.(^n)&^.,,(^n')e^""'e-^-"'^'e-^(-'"+'^™')(^-^') 

X Wm*(R - Xi)w™,(R' - Xj)G'a(n, k; iu;„)G'a(n', k'; ium') 

X Xn,k(R)X^„:k(R')Xn',k'(R)Xn',k'(R')- (25) 

The integrals over the imaginary time variables r and r' can now easily be performed, 
resulting in the conditions Un = ujn> = ujm + ^m'- Furthermore, if the optical lattice is 
deep, it can at every site be well approximated by an harmonic potential, characterized 
by the frequency uj. Then, we are also allowed to use the so-called tight-binding limit, 
in which the Wannier functions Wn are replaced by the tight-binding functions 0n, i-e. 

U7„(R-Xi) ~0„(R-x,). (26) 

Here, 0n(R) are the eigenstates of the three-dimensional isotropic harmonic oscillator 
potential given in spherical coordinates by 

-1/2 



^n(R) = 0n,£,m(R) 




+ 1/2 



n 



v/(^+l/2)! 



X e-^^/2'^4V2+^)((i?//)2) ^Rliy Yem{e,<P), (27) 

where L^^^'^^^\X) are the generalized Laguerre polynomials and Yi „ii(^, 0) the spherical 
harmonics. For the atoms, the harmonic oscillator length / is given by /a = \J'h/ma,uj] 
for the molecules, the harmonic oscillator length / is given by ls./\/2. Furthermore, in 
the extreme tight-binding limit the energies en,k in equation (PHj) are given by en,k = 
= (2n + 1 + 3/2)huj and only depend on the radial and angular momentum quantum 
numbers n and i, respectively. Similarly, we have that Ga(n, k; ium) = Ga(n; iujm)- Also 
note that in the tight-binding limit, the overlap for wave functions on different lattice 
sites is negligible. As a result, we can consider only the terms on the right-hand side 
of equation ()25p that have all atomic and molecular Wannier functions centered around 
the same site Xj. Furthermore, because of translational invariance of the lattice, we can 
set Xj to zero in the integrant. As a result, the right-hand side of equation (j^^ does not 
depend on k, k' anymore, and the sum over the lattice momenta yields the factor N^. 
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Next, we evaluate the sum over Matsubara frequencies, keeping the external 
frequency Um + ^m' = "^n of the molecules fixed. First, we rewrite the product of 
Green's functions as 

Ga(n, tiu„,)G^{n , tujrn') = -ri—— ■ — ■ — . (28) 

-in[uJm + CU^/) + en + en' - 2/i 

Subsequently, we perform the sum and take the two-body limit, which amounts to taking 
the atomic distribution functions, obtained after the Matsubara summation, equal to 
zero We find 

1 v.^ v.. V 



- -Tr [GaSGaS] = E E E E / d^d^ b*^^^Mk 



m,m' n,n' 

km* / 



-ihujn + en + en' - 2/i 
where 0n''^(-R') denotes the molecular /atomic tight-binding function. For systems of 
interest to us, the molecules occupy only the lowest Bloch band, i.e., we take the band 
indices m and m' of the molecular fields in the last expression equal to zero. The 
tight-binding function is given by 0}]^(R) = {2/7rlD'^^'^ e~^^^'"'^^ . 

Next, we need to evaluate the integrals over R and R' of the form 

J rfR0r(R)0n(R)0n'(R) = 2 ^ 

X [dRde d<i> sin^ ^-2R^L\^I'^'KR')lT^'\R') 

J 7(£+l/2)!y'(f + l/2)! 

xYt,m{0,(p)Ye,,„^,{d,(j))\ n I \ n' I ' 

where on the right-hand side of the above equation the integration variable R is made 
dimensionless using the length Z^- The integrations over the angles can be directly 
evaluated using the orthonormality relations for the spherical harmonics. To evaluate 
the remaining integral over R, we make use of the following relation 

'2\ 



r dX e-2^'x2+2^4^/2+^)(x2)Li!/'+')(X2 
Jo 



1 rdye-V''-'L^:^'^Hy)Lr^'\y) 

2 Jo 

T{n + n' + £ + 3/2) 1 



(31) 



2n\n'\ 2"+"'+^+3/2- 
Using this, we find that equation can be written as 

/ 2 T{n + n' + i + 3/2) 



J rfRC*(R)C(R)0?.'(R) = (-1)"^^- 



vr/? / n\n'\ 



a 



1 1 / n + £+l/2 \ n' + ^+l/2 \ 



- 2n+n'+£+3/2 (£+ 1/2)1 \^ n J \ u' J ' ^^^^ 

Finally, we have to perform the sum over the quantum numbers of equation ((221) • We 
first partially evaluate the resulting sum over the quantum numbers n and n' by summing 
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out all the contributions that have the same energy = e^^o = ^n/ + ^n'/ — Shu/2. This 
yields 



E 



rfR0r(R)0n(R)0n'(R) 



2 



6,,n+n'+i = 0r (O)C(O)' (33) 



where the relative wavefunctions of two particles in a harmonic potential ^^'^(r) will be 
introduced in section El and are given by equation (jHSj)- The reason for expressing the 
sum in terms of the relative wavefunctions 4'q'^{r) is to reveal the complete agreement 
between the field-theoretical calculation in this section and the two-body calculations 
in sectional 

For the second order term corresponding to the molecular self-energy we now obtain 
- ^Tr [G.SGaS] = -2g'Y.bhMKM 

i,n 

j^aa*/QWaa/Q\ 

^ ~ihujn + eq + ?>nuj/2-2^i ^ ^ 

This sum contains an ultraviolet divergence resulting from the use of pseudopotentials. 
In section El we show how to deal with this divergence by means of a renormalization 
procedure, which leads to the final form of the molecular self-energy /iSm 

, . 2^(i;iu;„-3/i^/2 + 2/i) 

nj:^{thujn) = g , (35) 

\'27il-lnuj 

where T{z) is the ratio of two gamma functions 
_ T{-z/2huj + 3/A) 
^^'^ = n-z/2hu + l/4y 

Note that the self-energy of equation (j35|l is completely equivalent to the self-energy 
obtained from a two-body calculation in ref. |13j . 

2.4- Generalized Hubbard model 

So far, we have calculated the self-energy of the bare molecules, but we have not yet 
taken into account the effects due to the purely molecular action and the purely 
atomic action S^- This is what we do next. 

The purely action S^, given by equation 0, can be rewritten by using equation 
which yields 

S,^= rfry rfx^;;(x,r) - + 5b -2/i + 2\/o(x)j ^/'^(x,r) 

= EE E Ci(^n)^m'j(^n) / C?> 
n i.i m.m' 



IX W^*(X - Xi 



i,j m,m' 

^2v72 



X (-ihio,, - ^ + 2K)(x) + <5b - 2/x ) <,(x - X,), (37) 



^a 



From which we can immediately read off the hopping term and the on-site energy 
em- For our purposes, the important case is when the molecules are in the lowest band. 



Dressed molecules in an optical lattice 
Then, we find for tlie flopping parameter 

Am. 



10 



tm = - J dx WQ*ix - Xi] 



+ 2Vn(x) 



Wn X — X 



wfiere i and j are nearest-neighbouring sites. The on-site energy em is given by 



dx. Wq^*(x — Xj 



4ma 



+ 2Voix) 



W'^ix - Xi). 



(3^ 



(39) 



As a resuit, the action Sm can be written as a iattice action 



(40) 



where {i^j) denotes the sum over nearest neighbours and em = ?iliu}/2 in the tight- 
binding fimit. 

fn exactfy the same way, we can rewrite the atomic action S^, for the iowest Bfoch 
band, yieiding 



with the atomic hopping parameter ta given by 



t^ = - j dx wl*{x - Xi 



2m„ 



+ Voix] 



wlix - X 



where i and j are nearest-neighbouring sites. For the on-site energy ea we find 



dx w^{x — Xi) 



2m„ 



+ K)(x) 



Wr,[X-Xi 



The hopping parameters ta.m satisfy exactfy the reiation 



t. 



huoX 



^(A/V27r«a,m)^ 



(41) 



(42) 



(43) 



(44) 



with A the waveiength of the fattice iaser and Zm = ^a/ "\/2- Furthermore, in the tight- 
binding fimit, ea = ^fiuj/2. 

Combining equations (jHll), (jHHjl and (finjl . we find that the on-site bare mofecufar 
propagator is given by 

— hG^ {ifiuJn) = —ifiujn + 'ihuj /2 — 2yU + 5 + hJ^miihiUn) ■ (45) 

with 6 the renormaiized detuning (see section Oj). We perform an anaiytic continuation 
ihujn E + iO and the zeros of the above equation are the poies of the Green's function. 
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which in turn correspond to the physical eigenstates of the on-site problem. To study 
these eigenstates, we calculate the spectral weight function p{E), given by 



(46) 



where G^\E) is the retarded Green's function of the system, i.e., G^\E) = G^{E+iO). 
In the two-body limit the chemical potential is zero and, because the self-energy is real, 
the spectral weight function becomes a set of delta functions located at the solutions e^- 
of the equation 



-hu + S = €0 



9 



T(e<, - 3huj/2) 



where we note that the delta functions have strength Z^-, given by 

\ -1 



1 



dhT.,^{E) 



dE 



(47) 



(4^ 



E = ta 



We now have all ingredients needed for the effective atom-molecule theory, since we 
can substitute G~^ and G^ and g in equation (jH)). In this manner, we have incorporated 
the on-site two-body physics exactly in the effective many-body theory. Furthermore, 
from the derived expression for the on-site bare molecular Green's function, we introduce 
the convenient notion of dressed molecular fields 60-, by making use of 

- {hh*) = GU^^co^) = E 7^ T- = -T. ZAbabl). (49) 

Using these dressed molecular fields, we finally find the Hamiltonian representation of 
our effective action, given by 

+ E E - 2/^) blA,<^ + E(ea - /i)4a, 



+ ^'EEV^. 



(50) 



denote the 
denote the 



where {i,j) denotes a sum over nearest-neighboring sites, aj/ctj 
creation/ annihilation operators of a single atom at site i, hl^^/h^^^ 
creation/annihilation operators corresponding to the dressed molecular fields at site 
and the effective atom-molecule coupling in the optical lattice is given by 



9 =9 



^o(x)0^(x)] 



(51) 



From equation (jSO)), we see that the effective Hamiltonian depends on several 
parameters. In particular, it depends on the energies of the dressed molecular fields 
and the wavefunction renormalization factors Z^- In this section, we calculated these 
parameters from first principles by a field-theoretical calculation. But from a comparison 
with ref. jTHj, we see that exactly the same results can be obtained from just a two- 
body calculation on a single site. This is as expected, since the poles of the molecular 
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Green's function correspond to the physical eigenstates of the on-site two-body Feshbach 
problem. 

So far, we did not take into account background atom-atom scattering, which is 
expected to be of particular importance for ^Li, which has an extremely large background 
scattering length. Because we use an effective Hamiltonian and dressed molecular fields, 
we can include the background scattering in a very convenient way. Namely, if we are 
able to solve the two-body physics on a single site including background scattering, then, 
through the notion of our dressed fields, the background scattering gets automatically 
incorporated into the many-body theory. This is what we achieve in the next section. 



3. Two-body physics on a single site. 

Consider the Feshbach problem for two atoms on a single site of an optical lattice. In 
the typical case of experimental interest, the optical lattice potential is deep and the 
energy of the atoms is low. As a result, the on-site potential is well approximated 
by an isotropic harmonic potential with angular frequency u. Furthermore, the atoms 
interact with each other through the potential Vaa which depends only on their relative 
coordinate r, while the atom- molecule coupling is given by V^m- The resulting problem 
can be separated into a center-of-mass part and a relative part. Only the relative part 
is important in solving the two-body physics on a single site, while the center-of-mass 
part determines the tunneling of molecules between adjacent sites. 

Taking all this into account, the relative Schrodinger equation for the two-channel 
Feshbach problem on a single site becomes 

I Km 5b a l^m) j [ l^m) ) ' ^ ^ 

where 6b denotes the energy of the bare molecular state and the relative noninteracting 
atomic Hamiltonian Hq is given by 

Ho = + -m^uj'r\ (53) 

ma 4 

Note that in writing down equation (j52p we have assumed that the energy of the bare 
molecular state \'4)m) is not affected by the optical lattice. This is well justified, since 
the spatial extent of the bare molecular wavefunction, centered around r = 0, is very 



small compared the the harmonic oscillator length Za = yh/ma^uo. Furthermore, we can 
rewrite equation ()52j) to obtain the following equation for the energy eigenvalues 

(^m|Km^^ 7^ T^KmlV'm) = E - 5^- (54) 

tj — Hq — Vaa 

3.1. No background atom-atom scattering. 

For many atoms of interest, such as for example rubidium or potassium, we have that 
|Ka| "C fiuj, which means that we can neglect Ka compared to Hq. The resulting 
problem was solved in ref. ^S] and here we will highlight the most important results. 
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The eigenf unctions of Hq can be written in terms of the generahzed Laguerre 
polynomials 

-r2/4«2 rl/2/ 2 /272N 

which correspond to the eigenenergies 

Er^ = (2n + 3/2)huj, (56) 

for n = 0, 1, 2, .... Using the completeness relation of the eigenfunctions \(f)^), equation 
(j54p can be rewritten as 

Invoking the pseudopotential approxiation for fermions, which means that (r| Vam|'?/'m) = 
g5{r), yields 

^(O)0-(O) 



E — En, 

m 

(5^ 



T(E) , m 
— lim 



where the function T(_E') was introduced in equation (j^ . The energy-independent 
divergence in equation ()58|) is an ultraviolet divergence resulting from the use of 
pseudopotentials, and was first obtained by Busch et al. in the context of a single- 
channel problem It can be dealt with by the following renormalization procedure. 
We define the renormalized detuning 6 as 

5 = (5b - lim-^, (59) 
'"^0 Ann r 

which has two major advantages. Not only do we absorb the ultraviolet divergence in the 
definition of S, but this renormalized detuning also has a relevant experimental meaning, 
in contrast to 6b- Namely, S corresponds to the detuning from the magnetic field Bq, at 
which the Feshbach resonance takes place in the absence of an optical lattice. This can 
be understood from the treatment of the homogeneous Feshbach problem without an 
optical lattice [21]. Here, the condition for the location of the Feshbach resonance is that 
the dressed molecular energy is equal to the threshold of the atomic continuum. This 
leads to the resonance condition 5b = limr^o^a5'^/47r^^r. As a result, the definition in 
equation (j59|) places the resonance in the absence of an optical lattice conveniently at 
5 = and, by construction, the (renormalized) detuning is of the following form 

5{B) = Afi{B-Bo), (60) 

where A/i is the difference in magnetic moment between the atoms in the open channel 
and the bare molecule in the closed channel. 

Substituting equation in equation (j3Hj) gives 

a^T(E) . , 
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which allows us to calculate the eigenenergies of the Feshbach problem on a single site 
given a certain detuning 6. Note the similarity between equation (jUT|) and equation (jlTj). 
The only two differences are a constant shift in the energy of 3hu/2 and a factor of 2, 
which are both readily explained. The shift of 3huj/2 is just the center-of-mass energy 
neglected in this relative calculation. The factor of 2 is due to the fact that here we 
consider fermions, while in the previous section we considered bosons. Since equations 
flTTj) and are for the rest completely equivalent, we see that indeed the energies e^- 
of the dressed molecular fields b^, needed for the effective Hamiltonian of equation (jHUjl . 
can just as well be determined from a two-body calculation, as from a field-theoretical 
calculation. An obvious advantage of the two-body approach is that it is much simpler. 
Furthermore, as we show in the next paragraph, the two-body approach also allows for 
an incorporation of atom-atom scattering effects, which is of course also possible, but 
much more difficult from a field-theoretical point of view. These effects are particularly 
important for the fermion ^Li, which has a very large background atom-atom scattering 
length Obg near the extremely broad Feshbach resonance at Bq = 834 G, namely on the 
order of — ISOOao, where is the Bohr radius. 

3.2. With background atom- atom interaction. 

The relative Schrodinger equation for two atoms in a harmonic potential interacting 
through the pseudopential = Vo6{r) without coupling to a molecular state 



can be solved analytically for s-wave scattering following a treatment along the lines 
of ref. j2Sl- We just note that in this reference a suspicious 'molecular' bound state 
appears, which we believe to be unphysical for several reasons. Most convincing is 
probably to consider the limit Obg in which case the interaction vanishes and we 

should recover the case of a simple harmonic oscillator. This is indeed what happens 
in ref. except for the suspicious bound state, whose energy goes to minus infinity. 
Since such a state does not arise in the treatment of the harmonic oscillator, we conclude 
that the state is unphysical and we will exclude it from further calculations. Note also 
that for the specific case of ^Li this issue is not relevant, since the suspicious state only 
occurs for positive values of a^g. 

The eigenvalues of equation (jU^ are given by E^, = (2v + 3/2)huj, which are the 
solutions to the equation 



where u is in general not an integer. The corresponding eigenf unctions (j)l^{r) = {r\(f)l^) 
have the following form 




(62) 



V2^ 



T{-E,/2huj + 3/4) ^ h_ 
T{-Ej2nuj+l/4:) ~ Obg' 



(63) 




(64) 



Dressed molecules in an optical lattice 



15 



where is a normalization constant and [/(— z/, 3/2, is a so-called confluent 

hypergeometric function of the second kind. We can determine A^, by applying a limiting 
procedure to theorem 7.622 of ref. [22], resulting in 

^2 ^ n-^-m .gg. 

" 4V2ttHIT{-u)[^o(^-u) -l/2)y ^ ^ 

in which denotes the digamma function. Furthermore, we can use the completeness 
relation of the eigenfunctions 10^*^) to rewrite equation which yields 

= E-6b. (66) 

In order to proceed from equation we realize that the integral (V'm|K,m|0^'^) 
only acquires a flnite contribution from a very small region around r = 0, because the 
spatial extent of both the molecular wavefunction and the atom-molecule coupling 
Vam is very small. This means that we are only interested in <i)^{r) for small values of 
r, giving 




27r/aA, + O - , (67) 



where in the second line we used equation ()63|) . Note that the eigenfunctions (\)'^{r) 
behave for small r completely analogous to the s-wave scattering states in the absence 
of an optical lattice ipkif ), which are of the following form 

. (68) 

kr 

where k is the relative momentum of the scattering state and 5o{k) is the so-called 
s-wave phase shift, given by 



k (69) 

5o(/c) = tan^"'^(— /cttbg). (70) 
Indeed, for small r, we have 

^,(r) = cos5o(A;) + + o{kr) 

rCV 

= -^^f^-l + 0(Pr)V (71) 
k \abg r J 

where in the second line we used equation ()70|) . As a result, for small r 

kip. A 



sm Oo{k) 



This relation is illustrated by flgure |21 where (pf,^, ip'f^ = y2Tikls,A^il)k/ sin^o and their 
difference are plotted as a function of r. We see that the difference indeed vanishes for 
small r and increases with increasing r. The scattering states ipkif) are the solutions of 
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r/4 r/4 

Figure 2. (a) Plot of the wavefunctions (pf,^ and -0^ as a function of the interatomic 
distance r. On the scale of the figure, the two wavefunctions cannot be distuinguished. 
(b) Plot of the difference between the wavefunctions (j)'^^ and ip'f. as a function of 
the interatomic distance r. Note that this difference is extremely small compared 
to the value of the wavefunctions themselves. All quantities on the axes are made 
dimensionless using the atomic harmonic oscillator length Za 



the relative two-atom s-wave scattering problem without optical lattice and without 
atom- molecule coupling [21]. The physical reason for the similarity between 0;y(r) 
and ipki'^) near r = is due to the following: for small r compared to /a? the atoms 
experience an effectively constant harmonic potential by which they are not affected. 
Therefore, for small r the wavefunction from the theory with optical lattice should, up 
to a normalization, reduce to the wavefunction from the theory without optical lattice. 
Since it is known how to solve the homogeneous Feshbach problem without an optical 
lattice 1211, ^t)le to profit from this knowledge. This is what we do in the next 

paragraph. 

3.3. Solving the two-atom Feshbach-problem. 

By combining equations and (ff2|) . we obtain 

^ 2nkXAl\{^P^\V^J^P,)\' 

l^sin^ 6o{k) E-E. =^-^^- 

From ref. [2lj, we have that 

(V-mlKml^fc) = (74) 

Substituting equations (jHSl) and (fTij) into equation (f7^ and rewriting the result with 
the use of equation (fTO)) yields 

? - {2u + 3/2)] [^o(_^) _ ^o(_^ _ 1/2)] = ^ - ^B. (75) 

Just like in section 3.1, we would like to renormalize the bare detuning ^b, which 
has no experimental meaning, to the renormalized detuning S, which gives the detuning 
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from the magnetic field Bq at which the Feshbach resonance takes place in the absence 
of an optical lattice. From the homogeneous theory that takes into account [23], we 
find at resonance that 

5b = (76) 
Ann |abg| 

Defining the renormalized detuning 5 as 5 = 5b — fi'^^a/47r/i^|abg| places the Feshbach 
resonance in the absence of an optical lattice conveniently at 5 = and brings equation 
(|75p into its final form 

2vrabg/^ ^[E- {2v + 3/2)] {i,%-v) - i,\-v - 1/2)] 

where by construction 6 = Afi{B — Bq) again. As in section 3.1, the term on the 
right-hand side of equation (f77j) can be interpreted as a molecular self-energy hT,^{E). 
Furthermore, upon going to a many-body theory, the solutions of equation (f77|) can be 
interpreted as the energies e„ of the dressed molecular fields b„, resulting in a generalized 
Hubbard Hamiltonian in the same way as before. 

For ^Li, which has a very large Obg, we expect the effects of the background 
scattering to be of particular importance. Therefore, in the next paragraph, we study 
equation (fTTjl for the specific case of ^Li to see if it leads to significant differences from the 
results given by equation (jHH), which doesn't take background scattering into account. 



3.4- Lithium 

Figure shows the energy-levels of two interacting ^Li atoms in the hyperfine states 
|1) and |2) on a single site near the Feshbach resonance at 834 G. The diagram has been 
obtained by numerically solving equation (|77|). where we used the known experimental 
values for g{B) and abg(-B) corresponding to ^Li Note that for the extremely 

broad Feshbach resonance at 834 G these parameters depend on the magnetic field B. 
Furthermore, we took for /a the realistic value 10~^ m, corresponding to u = 10^ s~^. In 
figure the same diagram is shown, but this time without taking the background atom- 
atom scattering into account. This diagram has been obtained by numerically solving 
equation ()61|) . where again the experimental value of the parameter g{B), corresponding 
to ^Li, was substituted. 

From figure El we see that for ^Li near the Feshbach resonance at 834 G it is 
important to include background atom-atom interactions, since it leads to a significant 
adjustment of the corresponding two-body energy-level diagram. Upon inclusion of the 
background atom-atom interactions, the horizontal asymptotes are shifted to different 
energy values and the avoided crossings, see also ref. [THj, become much less broad. 
Physically, this is a result of the fact that the large value of Obg effectively reduces the 
atom-molecule couphng. 
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Figure 3. (a) Relative energy levels (solid lines) as a function of the magnetic 
field B for a system consisting of two ®Li atoms in the hyperfine states |1) and |2) 
near the Feshbach resonance at 834 G in a harmonic potential oi uj — 10^ s~^ with 
background atom-atom scattering taken into account. The dashed line corresponds 
to the binding energy of a dressed ^Li2-molecule in the absence of an optical lattice. 
Note that as the eigenenergy of the dressed molecular ground state with optical lattice 
(solid line) decreases, it converges to the dressed molecular eigenenergy without optical 
lattice (dashed line), (b) Similar to panel (a), only this time the background atom- 
atom scattering is not taken into account. Note that as a result both panels differ 
significantly. 

The difference between equations (jUT|) and ()77p in the case of ^Li becomes even 
clearer when we look at the corresponding wavefunction renormalization factors given 
by equation ()48|). which are needed for the construction of the effective Hamiltonian 
given by equation (jHUjl . For the calculation of Z„ we need the molecular self-energy 
^Sm, which is given by the right-hand side of equation ()6ip for the case without 
background atom-atom scattering and by the right-hand side of equation (|77jl for the 
case with background atom-atom scattering. In figureEJ the corresponding wavefunction 
renormalization factors Z^^ are plotted for both function of the magnetic 

field B. Here, Zq corresponds to the groundstate of the relative on-site two-body 
Hamiltonian, Zi corresponds to the first excited state, and so on. Note that the two- 
body renormalization factors Z^r obtained including background interactions (figure^), 
are very different from the factors obtained without background interactions (figure EJd). 
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Figure 4. (a),(b) Wavefunction renormalization factors Zcr for ^Li near the Feshbach 
resonance at 834 G as a function of the magnetic field B with background atom-atom 
scattering taken into account. Zq is the renormahzation factor corresponding to the 
ground state of the two-body on-site Feshbach problem, Zi is the renormalization 
factor corresponding to the first excited state, and so on. (c),(d) Wavefunction 
renormalization factors Z„ for ^Li as a function of the magnetic field B without 
background atom-atom scattering. Note that the panels on the left differ significantly 
from the panels on the right. 

Also note that these factors are directly experimentally observable as demonstrated by 
Partridge et al. [7]. We thus conclude that background atom-atom scattering plays a 
significant role for an atomic gas of ^Li atoms in an optical lattice near the Feshbach 
resonance at 834 G. 

4. Conclusion 

In summary, we have derived from first principles an effective atom-molecule theory 
that describes an atomic gas in an optical lattice near a Feshbach resonance. The 
theory was formulated in the tight-binding limit and incorporated all the relevant two- 
body physics exactly. The field-theoretical derivation reconfirmed that an atomic gas 
in an optical lattice near a Feshbach resonance is accurately described by a generalized 
Hubbard model, as first obtained in ref. ^Hl- In the original approach of ref. 112] the 
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generalized Hubbard model was derived in an easier way by starting from a two-body 
calculation on a single site. Due to the equivalence of both approaches, we used the 
latter to formulate a way in which to include the background atom-atom scattering into 
the many-body theory, which has not been done previously. To this end, we showed 
how to solve exactly the two-channel Feshbach problem for two atoms on a single site 
including background atom-atom scattering. The solution was applied to ^Li near the 
experimentally relevant Feshbach resonance at 834 G. Specifically, the two-body on-site 
energy levels and the wavefunction renormalization factors Z^r were obtained, which are 
needed for the generalized Hubbard Hamiltonian describing the many-body physics. As 
it turned out, the various energy levels and the wavefunction renormalization factors 
obtained with background interactions taken into account, are significantly different 
from the energy levels and the renormalization factors obtained without background 
interactions taken into account. From this we conclude that in the case of ^Li background 
atom-atom scattering plays an important role and cannot be neglected for an accurate 
microscopic description of future experiments with atomic lithium in an optical lattice. 
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